Brief remarks on "Electric dipole moment enhancement factor of thallium' 
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In a recent paper by Porsev et al. [arXiv:1201.5615vl], the authors have claimed to have resolved 
the controversy arising from the different ab initio results available for the EDM enhancement factor 
of Tl. In our opinion, any such attempt to resolve the discrepancies between different calculations 
has to thoroughly compare the basis sets used, methods employed and approximations considered 
in the two different cases. However, Porsev et al. have not succeeded in doing so in their current 
paper. We clarify some of their misunderstandings about our work and address some specific issues 
in this note. 



First of all, we would like to state that some of the 
points made by Porsev et al. in their paper [l[ are factu- 
ally incorrect. We quote, for example, their statement: 
"Since the valence- valence correlations are very large, the 
CI method provides better description of these correla- 
tions than the perturbative approaches such as RCC due 
to possible large contributions of higher-order (or higher- 
excitation) correlations.". The fact is that the coupled- 
cluster (CC) theorj^or its relativistic extension, RCC, is 
NOT perturbative 0,11 even though it is equivalent to all 
order perturbation theory. Furthermore, the configura- 
tion interaction (CI) approach is a subset of CC method 
at the same level of hole-particle excitations and hence, 
it does not include any higher-order correlations beyond 
CC theory as they have claimed. For example, CC sin- 
gles and doubles (CCSD) contains all the terms that are 
present in CI singles and doubles (CISD) and much more 
[3,3. Further, we would like to clarify that the valence- 
valence correlations are very large in thallium (Tl), as 
they have reported, only if one treats 6s^ 6p as valence 
electrons. However, that does not preclude one from con- 
sidering 6s^ as part of the core and 6pi/2 as valence and 
treating the valence-core interactions to all order in the 
residual Coulomb interaction via the CC theory as we 
have done Therefore the valence- valence correlation 
of Porsev et al. [H and Dzuba and Flambaum (D & F) 
is a part of our valence-core correlation. 

Furthermore, the same many-body approach is em- 
ployed by Porsev et al. [l| for their calculations using 
V^~^ as well as V^~^ orbitals. In both cases, they 
have treated the three outermost electrons as part of the 
valence space using CI and the rest of the electrons as 
core using MBPT. In this context, they mention that the 
RCC approach employed by us does not treat 6s 6p^ and 
6s^ ns, where n being > 7, on an equal footing. However, 
this is not true. We have treated both of them on the 
same footing, that is, as opposite parity excitations from 
a common reference state [4|. 

Porsev et al. [l| are correct in stating that: "It is very 
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important to accurately account for the contributions of 
the 6s 6p^ configurations." . However due to their incor- 
rect assumptions that we have mentioned above, they 
have arrived at an incorrect conclusion that "In the RCC 
method (in referring to our work in [3]), these contri- 
butions were treated as excitations of the core electrons 
which is unlikely to provide the required accuracy." They 
imply that 6s^ electrons has to be treated only as valence, 
or in other words, one has to treat Tl only as a three- 
valence atom as against our treatment of it as a mono- 
valent atom. Referring to Tl as a one-valence or three- 
valence system is a matter of semantics. However, what 
really matters is not the terminology one uses (three- 
valence approximation vs. one- valence approximation), 
but rather the physical effects contained in a particular 
theory. Indeed, Porsev et al. [l| and D & F [5| have 
considered 6s^ and 6pi/2 as valence electrons and have 
evaluated the valence- valence correlations by CI. The net 
level of excitation of the configurations in their CI cal- 
culations does not appear to be beyond triples with ref- 
erence to 6s^6pi/2- In contrast, our CC theory has all 
linear and nonlinear single and double excitations from 
all the core and valence electrons (in particular 6s^ and 
6pi/2 electron) and their products as given in Eqs. (3 & 
4) of our paper Q. In other words, in addition to sin- 
gle and double excitations from the 6s^6pi/2 Fermi vac- 
uum state, our wave function consists of a large number 
of triple, quadruple, quintuple, sextuple, septuple and 
octuple excitations which are obtained as disconnected 
products of single and double excitations. Therefore in 
addition to all the configurations included in the calcu- 
lations of Porsev et al. [l[ and D & F [Hj , our calculation 
contains many more configurations corresponding mostly 
to higher order excitations. From the information given 
on the configurations included in the model space and 
the excitations considered for the CI calculation of D & 
F, it does appear that this calculation is not size exten- 
sive [3, S] as the CI is not complete in the chosen orbital 
space. Even though Porsev et al. do not specify their 
configurations in the model space, it seems as though 
they have not performed a full CI, particularly in view 
of their statement "higher n orbitals were allowed fewer 
number of excitations" , and therefore it also might not 
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be size extensive. In contrast, our RCC calculation is size 
extensive at all levels of excitation. The probable lack of 
size extensivity in the two above mentioned CI calcula- 
tions could be one of the likely sources of discrepancy 
between those calculations and ours. 

Further, we also would like to clarify that contrary to 
the opinion expressed by Porsev et ai, the size of our 
basis set is not small. If our basis functions were not 
nearly complete then we would not have obtained good 
results for both the electric dipole (El) amplitude and 
the magnetic dipole hyperfine constants shown together 
in Table IV of our paper Very surprisingly, Porsev 
et al. make no reference at all to these results. The 
wave functions at both near- and far-nuclear regions are 
subjected to scrutiny through these properties. On the 
whole our results are in better agreement with experi- 
ments than those of D & F's and Porsev et al.^s results; 
For example, our RCC method gives El amplitude of 
the 75 — > 6P1/2 transition as 1.82 au against the experi- 
mental value 1.81(2) au. In contrast, this is reported as 
1.781 au and 1.73 au by Porsev et al. and D & F 
0, respectively, which are outside of experimental limit. 
The good agreement with experiments of our results for 
the above mentioned two properties which are related to 
the EDM enhancement factor calculation is a strong in- 
dication that our EDM enhancement factor is reliable. 
In this context, we would also like to mention that we 
have considered 38s; 34pi/2;3/2; 34d3/2;5/2; 30/5/2;7/2 and 
20(77/2;9/2 number of orbitals for the SCF calculation. 
The number of basis functions and the basis parameters 
in the finite basis space are well optimised by comparing 
the single particle energies and SCF energy with respect 
to the numerical (non-parametrical) results obtained by 
the general-purpose relativistic atomic structure program 
(GRASP). We have performed various tests at the Dirac- 
Fock level to check the correctness of our single particle 
orbitals like the comparison of single particle energy dif- 
ferences, dipole and EDM matrix elements with those 
obtained using GRASP (details can be found in Q). Af- 
ter the SCF step, we have truncated the virtual space 
by dropping the high-lying virtuals whose contribution 
to the results is small. As a matter of fact, those basis 
fimctions which are used by Porsev et al, B-splines used 
by D & F and Gaussians used in our work have different 
qualitative behaviours at different regions. Therefore, 
one cannot directly compare the number of basis func- 
tions used in different calculations. Having said that, we 
would like to emphasize that the use of Gaussian basis 
functions in relativistic calculations have been well tested 
for a number of different properties. 

We also emphasize that the EDM enhancement fac- 
tor for Tl depends mainly on the following three factors: 
the EDM matrix elements, the El matrix elements and 
the energy differences between the ground state and the 
intermediate states. Porsev et al. have calculated the 
6pi/2 ~ 75 EDM matrix element using the CC method 
restricting the number of basis to n = 14 for all partial 
waves in an attempt to reproduce our result and have re- 



ported that it reduces the value by 18% when compared 
to their calculation with a bigger basis set considered for 
the rest of the calculations. However, this reduction will 
not be reflected in the total EDM enhancement factor 
by a similar magnitude. As the energy differences be- 
tween the high-lying virtual states and the ground state 
appearing in the denominator will be large, the contribu- 
tion from those states to the overall EDM enhancement 
factor would be small. 

Porsev et al. have calculated the core-valence correla- 
tions using the CI-|-all-order method only in the V'^~^ 
case and they report that their contribution to the total 
EDM enhancement factor is less than 1%. Further, they 
assume that the magnitude of these correlations will more 
or less be the same in the case. As the nature of 

the orbitals are very different in the two cases, for details 
please see 0], this assumption may not be valid. In Ref. 
[7| we have demonstrated by numerical calculations that 
the EDM enhancement factor of Tl computed using the 
V^-^ orbitals over-estimates the result in comparison 
to that calculated using the V^~^ orbitals at the Dirac- 
Fock level of the theory. Further, their CI-|-all-order ap- 
proach only includes the linearized CC terms and hence, 
the contributions from the non-linear CC terms which are 
omitted in their calculations may also have non-negligible 
contributions. 

Porsev et al. have used the sum-over-states approach 
to calculate the EDM enhancement factor from the 
two specific states: 6s^ 7s and 6s^ 8s using the RCC 
method and they compare those results with their full 
Cl-f MBPT-f RPA results and both these results agree to 
each other, within 2% accuracy. Further they have com- 
mented that the result inferred from the Figure 2 of our 
paper is 10% lower than their RCC result for the 6s^ 7s 
state. However, we would like to remind the fact that we 
had given the combined contributions of the following 
RCC terms: DT^^\ DS[l^, s[°J''DS[1^ and sil^''DS[l^ 
to the singly-excited intermediate states of s symmetry 
in Figure 2 of our paper Q , as quoted in the text therein. 
It is not clear to us whether or not Porsev et al. have 
included the contributions from all these RCC terms in 
their results and if not, then the comparison may not be 
very meaningful. 

On a note aside, we remark that the calculation of the 
effective electric field of a molecule such as YbF is differ- 
ent from the calculation of the EDM enhancement factor 
of Tl on the following grounds: the question of one- or 
three-valence does not arise in the former case and in 
addition, the amount of correlation to the effective elec- 
tric field of YbF is very small, about 3%, as observed in 
all the earlier calculations including our own recent RCC 
calculation, unlike the case of Tl where the correlation 
contributions are large. So, we do not agree with the 
view expressed by Porsev et al. that the YbF calcula- 
tions are more difficult than the Tl calculations. 

In summary, the current attempt by Porsev et al. [l| 
in trying to resolve the discrepancies between different 
EDM enhancement factor values for Tl has merely added 
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a new result to the literature. Although, it reports an 
agreement, within the limits of the quoted uncertain- 
ties, with the earlier two calculations 0, il, it clearly 
emphasizes that the correlation effects ignored by D & F 
in their CI+MBPT+RPA calculation could be very im- 
portant and it further quantifies the magnitude of these 
corrections to be as large as 7% in the V^~^ potential 
itself, despite the severe accidental cancellations among 
the different corrections. Porsev et al. also observe that 
the energies and various properties computed using the 
orbitals generated in the V^~^ potential agree quite well 
with the experiments when compared to those results of 
y7V-3 potential in concurrence with our observation re- 
ported in [7]. Fortuitously the net result of Porsev et 



al. agrees with the results of fs', IsJ . We observe that the 
method followed in the work of Porsev et al. [l| is not 
very different from that of D & F Q and their RCC calcu- 
lation of the EDM enhancement factor is not as complete 
as ours. In the absence of adequate information about 
various calculations performed by Porsev et al., particu- 
larly their CI calculation and a detailed comparison be- 
tween their intermediate results and ours given in Table 
I of our paper Q], it is difficult at this stage for us to 
judge their work. However, due to the reasons discussed 
in this note we believe that Porsev et afs work on the 
EDM enhancement factor calculation of Tl as reported in 
m has come nowhere close to resolving the discrepancies 
between the different calculations. 
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